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ABSTRACT 

A consensus is emerging that interacting galaxies show depressed nuclear gas metallicities compared 
to isolated star-forming galaxies. Simulations suggest that this nuclear underabundance is caused by 
interaction-induced inflow of metal-poor gas, and that this inflow concurrently flattens the radial 
metallicity gradients in strongly interacting galaxies. We present metallicities of over 300 H II regions 
in a sample of 16 spirals that are members of strongly interacting galaxy pairs with mass ratio near 
unity. The deprojected radial gradients in these galaxies are about half of those in a control sample 
of isolated, late-type spirals. Detailed comparison of the gradients with simulations show remarkable 
agreement in gradient distributions, the relationship between gradients and nuclear underabundances, 
and the shape of profile deviations from a straight line. Taken together, this evidence conclusively 
demonstrates that strongly interacting galaxies at the present day undergo nuclear metal dilution due 
to gas inflow, as well as significant flattening of their gas-phase metallicity gradients, and that current 
simulations can robustly reproduce this behavior at a statistical level. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: interactions — galaxies: 
ISM 



1. INTRODUCTION 

The production and redistribution of heavy (non- 
hydrogen or -helium) elements within galaxies is an im- 
portant aspect of galaxy evolution. The amount of heavy 
elements in the gas phase of galaxies has been a topic of 
study over the decades since resolved galaxy spectroscopy 
became commonplace. In particular, emission line spec- 
troscopy of star-forming regions has been of great use 
in understanding the distribution of chemical elements 
in nearby, normal galaxies. For more distant galaxies, 
the mass- metallicity relationship, an important signpost 
of galaxy chemical evolution, has been used to constrain 
chemical evolution across a broader range of galaxy types 
and, increasingly, across cosmic history. 

Intense study of the mass-metallicity (hereafter, M — 
Z) re lation in recent years has yet to pin down i t s origin 
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How- 
ever, it is apparent from an increasing number of unique 
studies that part of the scatter in the relation is due 
to a nucle ar underabundance in interacting galaxie s, at 
both low ( |Lee et al.|[2004| |Ekta fc Chengalu r ]| 201 0D and 
high masses jKewley et al 2006a; Ru ple°CT°lu f ! J2008 



Ellison et "aL]|2008[ |Michel-L>ansac et al.||2008j |HeepIes| 



et al.||2009[ |Sol Alonso et al.||2010[ ). If one of the pro 

gemtor galaxies m an interacting syste m is a spiral, a 



radia l abundance gradient exists (e.g., Zaritsky et al. 
19941 Ivan Zee et al. 1998), with lower abundances at 



larger galactocentric radius. During the interaction, low- 
metallicity gas from the galaxy outs kirts is torqued into 
the high-metallicity galaxy center (Mihos & Hernquist 
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1996 Barnes & Hernquist 1996), resulting in gas with 



a lower" average abundance (as first suggested in Kew- 
ley et al. 2006a). We discussed the quantitative details 
ol this scenario for present-day equal-mass mergers in 
|Rupke et al.| ( [2010| . 

A prediction ot this model is that this gas redistri- 
bution dramatically flattens the initial rad ial metallicity 



gradient very shortly afte r first pericenter ( Kewley et al. 



2006a||Rupke et al.||2010|). |Chien et al.| |2007|) found no 



oxygen gradient along the tidal tails of the Mice, consis- 
tent with strong radial mixing of gas. Shallow or bro- 
ken gradients are also present in many nearby barred 
galaxies, which are undergoing radi al redistribution of 
gas due to bar-induced gas motions flVila-Costas fc Ed 
munds|[T992l |Zaritsky et al.|[L994) |Martin fc Roy||199 
Roy fc Walsh|ll997[ |Dutil fc Roy|ll999p . Barred galax 
ies also show lower central metallicity tha n galaxies of 



or 

% 



simila r morphological type and luminosity (Dutil & Roy 
19991, as is observed in mergers. Numerical models ot 



barre d galaxies reproduce this behavior (Friedli et al. 



1994) 



ma companion paper ( Kewley et al.|20l0 ), we showed 



that eight interacting galaxies which are underabundant 
with respect to the luminosity-metallicity relation all 
have gradients that are flatter than the nearby spirals 
M83, M101, and the Milky Way. This result strongly 
supports the model of merger-induced metal mixing. In 
the present paper, we expand on t his result by giving a 
detailed analysis of the systems in jKewley et al. ( 2010[ ), 
doubling the sample of interacting galaxies with mea- 
sured gradients (to 16), and presenting a comprehensive 
comparison to isolated spirals and the numerical simu- 
lations. The interacting galaxies in this sample are all 
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in the early stages of merging. They are after first pas- 
sage of the two (or more) galaxies, but probably prior to 
second passage. A future paper will present similar data 
on galaxies in the later stages of merging (Rupke et al. 
2010, in prep.). 

Understanding chemical evolution during a major 
merger will shed light on galaxy evolution at the peak 
of cosmic star formation, since the major merger rate, 
espec ially involving g^as-rich galaxies, was high er in the 
past ( |de Ravel et al.|[2009[ |Bundy et al.||2009[ ). These 
major mergers have likely played an important role in 
the producti on and growth of massive, red galaxies sinc e 



mag; Kochanek et al. 2001|), and have to tal infrared lu- 



minosities near L" IR ( Sanders et al.|2003 ). As mentioned 



1 (e.g., Bundy et al. 2009 Robaina et al. 2010) 



Simulations suggest that dry mergers (with n o gas) ha' 



ten stellar metallicity gradients in ellipticals (Di Mattco 
et al.|2009 ). It may be that wet mergers that form bulge 



dominated systems may do the same, if significant star 
formation occurs after the gas gradient flattens. In short, 
the metallicity redistribution caused by gas-rich mergers 
may significantly impact galaxy chemical evolution. 

In §§ 2 and 3, we present the observations and our anal- 
ysis. In § 4 we present the metal distributions, and then 
compare to numerical simulations in § 5. We summa- 
rize in § 6. Cosmological quantities are computed using 
H = 73 km s _1 , il m — 0.27, and ft = 1, and in the 
rest frame defined by the cosmic microwave background 
(where applicable). 

2. SAMPLE AND OBSERVATIONS 
2.1. Sample 

Our sample consists of galaxy pairs undergoing strong 
interactions in the local univer se. The pairs are cho- 
sen fr om both optically-selected (Arp 1966; Barton et al. 
2000| and i nfrared-selected (Sanders et al. 2003; Surace| 
et al. 2004) catalogs of interacting systems! The pri- 
mary selection criteria, to ensure that the systems are 
strongly interacting, are that the projected nuclear sep- 
aration is small (< 30 kpc), the galaxy mass ratios are 
small (1:1—1:3), and that the systems are obviously dis- 
turbed (i.e., after first passage). To distinguish from sys- 
tems in later stages of merging, we have rejected systems 
with (1) projected pair separations <15 kpc and/or (2) 
compact morphology and faint tidal features that indi- 
cate the system is near or after nuclear coalescence. The 
first criterion may remove some systems very near first or 
second pericenter, but also securely rules out most sys- 
tems after second pericenter, since rapid orbital decay 
in strong, equal-mass interactions leads to small nuclear 



separations after second pericenter (e.g., |Barnes|[l998 
and references therein) 

We have chosen spiral-spiral interactions to maximize 
the prevalence of H II regions. In each system in our 
sample, there is an obviously interacting pair with small 
projected separation, near-equal masses, and signs of in- 
teraction. However, in some cases there is also a third 
galaxy at large projected separations (>50 kpc) from the 
dominant pair which may or may not be involved in the 
interaction. Given their large distances, we do not fur- 
ther consider these possible third members of the system. 

Basic properties of individual galaxies in our sample 
are listed in Table [l] Our sample of 9 pairs/groups have 
fairly average properties for nearby galaxies (Figure fl]) . 
Individual galaxies are evenly distributed around L* K (±2 



above, the pair mass ratios (as traced by X-band lumi 
nosity) are in the 1:1—1:3 range, with projected nuclear 
separations of 15 — 30 kpc (except for Arp 248, where the 
separation of the two primary galaxies is 54.3 kpc). 

Table [T] also lists optical diameters, inclinations, and 
line-of-nodes position angles for the sample. Man y of 
these parameters are taken directly from HyperLeda ( Pa- 
turel et al.|[2003| ). For the NGC 2207 / IC 2163 system. 
we used parameters det e rmine d from numerical models 



( Elmegreen et al.| 1995 12000). In the cases of a few 
galaxies, the tidal stretching of the galaxies or general 
interaction-induced disturbances make determining disk 
orientations a challenge. In this handful of cases, we 
have turned where possible to H I maps to determine 
the position angles. These cases are listed in the table 
notes. We have also used this H I data and inspection 
of visual images to double-check HyperLeda inclinations. 
The inclinations of two galaxies remain uncertain: Arp 
248 NED02 and Arp 256 NED02, which are particularly 
tidally-stretched. In both of these cases, the effect on the 
gradient is compensated for somewhat by a larger mea- 
sured R25- We estimate that the inclinations lie in the 
range — 45°, and compute the deprojection using 30°. 

2.2. Control 

A control sample of isolated galaxies is a necessary 
component of this study. In particular, we need to com- 
pare our results to a group of galaxies that plausibly 
represent the galaxies in our sample as they were be- 
fore a major interaction. We have selected our control 
sample from published data on nearby, isolated spirals 
that have data on H II region emission-line fluxes. In 
particular, we have focused on mid-to-late-type spirals 



with a mix of bar strengths ( Martin 
1999), since mid-to-late types are bes 



1995] |Dutil & Roy 

represented in the 
literature. Gradients correlate with bar strength in a re - 
lationship with significant scatter ( Dutil & Roy |1999 ^ 
and depend we akly or not at all on spiral type (|Zarit 



sky et al. 1994). Our sample should thus be represen- 
tative ol isolated spirals at the level allowed by small 
number statistics. Finally, with only two exceptions, we 
have picked galaxies from moderate-size, homogeneous 
studies to minimize systematic uncertainties. (The two 
exceptions, NGC 300 and M101, are well-characterized 
benchmarks in gradient studies, and so we include them 
as well.) The salient properties of our control sample, 
and the references for the sample properties, are listed in 
Table H 

The interacting and control samples are well-matched 
in mass, as traced by -ftT-band luminosity (Fig. [Tji; the 
control sample is only half a magnitude fainter on aver- 
age). They are also well-matched in optical radius, with 
average (median) optical radii of R25 = 13 (12) kpc in 
both samples. However, the interacting galaxies have an 
average (median) infrared luminosity that is boosted a 
factor of 7 (4) above the control sample (Fig. II|)) . This 
is unsurprising, given the expec ted higher star form ation 
rates in the interacting sample ( Barton et al.|[2000 ). 

Rather than use gradient determinations published in 
the literature, we have used the published line fluxes 
to measure abundances using the same strong-line di- 
agnostic that we apply to our interacting galaxies sam- 
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pie Q3.2). We have also performed our own spatial de- 
projections using published photometric parameters (see 
references in Table p| . The resulting oxygen abundance 
gradients, in terms of dex/kpc and dex/R,25, are given in 
Table [2 

2.3. Observations 
We used the multiobject spectroscopic capabilities of 



Oke et al. 


1995 


McCarthy et al. 


1998 


) to observe H II 
be excellent blue 


regions in the galaxies ot interest. T 



sensitivity of LRIS was important for measuring the 
[O II]AA3727, 3729 emission lines at high S/N; the [O II] 
doublet is used in our strong-line metallicity diagnostics 
(jj3]). Our typical setup was to observe the wavelength 
range 3500 — 7000 A, using the 900 1/mm grating in the 
red and either the 400 or the 600 1/mm grism in the 
blue, and the 560 nm dichroic to separate the red and 
blue spectra. (For one galaxy, Arp 248, we were forced 
to use the 300 1/mm grism in the blue to obtain the 
entire spectrum, because LRIS-R was inoperable.) For 
over half of our data, the LRIS atmospheric dispersion 
corrector (ADC) was in place, mitigating any light loss 
in the blue. For the other half (taken prior to mid-2007), 
slit losses due to atmospheric refraction were generally 
small. Total exposure times per mask were typically an 
hour. 

Slitmasks were created in most cases using sensitive 
Ha emission-line images. Slits were 1" wide. The Ha 
images were taken using redshifted narrowband filters 
at either the Vatican Advanced Technology Telescope 
(VATT; E. Barton and R. Jansen 2005, private comm.) 
or the University of Hawaii 88-inch telescope. However, 
in two cases (Arp 256 and Arp 298), the spectroscopic 
data was taken in the context of a program to observe 
star clusters (Chien et al. 2010, in prep.). In this case, 
star clusters were targeted using Hubble Space Telescope 
(HST) broadband F435W images. Since significant Ha 
emission fell in the slits in these cases, we were able to 
extract useful H II region spectra for the current work. 

3. DATA REDUCTION AND ANALYSIS 

3.1. Data Reduction 

The LRIS data were reduced using a pipeline we devel- 
oped specifically for the current dataset. This pipeline 
combines Image Reduction and Analysis Facility (IRAF) 
and Interactive Data Language (IDL) scripts. For each 
side of the spectrograph, the steps we took were as fol- 
lows: bias subtraction; slit identification and tracing us- 
ing Hatfield exposures; arc lamp identification; slit ex- 
traction and spatial rectification (including any atmo- 
spheric dispersion correction, if required); a small shift to 
match dispersion solutions among exposures; data com- 
bination, including noise-threshold cosmic ray rejection; 
wavelength calibration; extraction of spectra; sky sub- 
traction; and flux calibration. In the slit identification 
step, slit edges were detected using a directional filter 
that computed the first derivative of the flatfield in the 
spatial direction. We extracted each local peak in Ha 
using a l"xl" aperture, and subtracted the sky using 
off-object spectra free of galaxy continuum and line emis- 
sion. Sky subtraction was generally imperfect due to the 
change of the LRIS point spread function across the 6' x 8' 



field-of-view combined with the use of a small number of 
sky spectra that were typically on the field edges (to 
avoid galaxy contamination within the on-source slits). 

Blue and red spectra were reduced separately, and then 
combined near 5600 A by multiplicatively matching the 
continuum fluxes in this region. Following the analysis 
of the spectra, this relative flux calibration was further 
refined. 

Spectra were mapped to positions on the sky assuming 
correct alignment of each slitmask during observations 
along with knowledge of the slit plane to focal plane 
mapping. (We used at least four alignment stars per 
slitmask to ensure correct mask alignment.) Galactocen- 
tric radii were then computed using info rma tion on the 
3-dimensional alignment of each galaxy ({2.1). To deter- 
mine the center of all galaxies except NGc 3994/5, we 
used near-infrared positions from 2MASS. Sloan Digital 
Sky Survey (SDSS) centers were used for NGC 3994/5. 
Though these systems are gravitationally disturbed, the 
near-infrared is able to probe the central bulge in each 
galaxy and thus yield a reasonably accurate true gravi- 
tational center. 

3.2. Spectral Analysis 

We extracted information from the resulting spectra 
by performing detailed fits to the stellar continuum and 
line emission. The software fitting package that we wrote 
for this purpose, UHSPECFIT, was based on th e stellar 
fitting routines of Moustakas fc Kennicutt| ( |2006[ ), who fit 
stellar continua using a linear combination of several tem- 
plate spectra. We built on these routines a sophisticated 
suite of emission-line fitting software that fits all emis- 
sion lines simultaneously and allows for multiple velocity 
components , start ing with the code framework used in 
Zahid et al-lpoTu] ). 

'lb fit these relatively high-resolutio n data, we used the 
stellar popul ation synthesis models of Gonzal ez Delgado] 
et al. (2005) as our stellar continuum templates. We 
relied on the solar metallicity, Geneva isochrone models 
for our fits. The quality of the fits is not a sensitive 
function of stellar metallicity; since our goal was not to 
constrain the stellar metallicity, but rather to remove 
the stellar continuum, we used only the solar metallicity 
models. The stellar fits were performed on the blue half 
of the spectra (from 3700 — 5500 A in the rest frame) 
after removal of wavelength regions near strong emission 
lines; the resulting stellar fit was then normalized to fit 
the red half of the spectra. We subtracted the stellar 
continuum from the spectra and performed emission-line 
fits. Two example fits are shown in Figure [2j One is 
a high emission line flux case, and the other is a low 
emission line flux case (roughly 1%, in Ha flux, of the 
high flux case). These are illustrative of the typical high 
fidelity of the fits, and reveal that fits to the Balmcr 
absorption lines in particular are very good. 

Once the fits were complete, a final blue-red relative 
flux calibration was performed by comparing the mea- 
sured reddening from the H/3 /H7 emission line ratio 



(using the extinction law of Cardelli et al.||1989[ ) with 
that measured from the Ha/H/i ratio. The discrepancy 
was used to correct the blue emission-line fluxes. This 
method worked very well for our data (which showed high 
signal-to-noise in H7) and our choice of spectral fitting 
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(which permitted robust separation of Balmer absorption 
and emission). However, our abundance errors increase 
with the use of H/3/H7 to determine the reddening (as 
opposed to Ha/H/3). 

We corrected the emission line fluxes for reddening 
( Cardelli et al.|1989 ) and computed important line ratios. 
Spectra with very uncertain reddening because of poorly- 
fit or weak H/3 and H7 emission lines were ignored. The 
dereddened line ratios were then used to classify each 
spectrum as H II-regi on, composite, Seyfert, or LINER 
(Kewley et al. 2006b). We calculated abundances only 
for those regions which are securely of H II-region ori- 
gin, without probable contamination from other sources 
of excitation. We included only those data identified 
as H II-region spectra based on the [N II]A6583/Ha vs. 
[O III1A5007/H/3 diagnosti c (i.e., regions that fall below 
the |Kauffmann et al |2003 star formation locus) and H II 
or composite spectra based on the [O I]A6300/Ha vs. 
[O III1A5007/H/3 diagn ostic (i.e., regions that fall below 
the Kewley et al.|2001 extreme starburst line) . For both 
line ratio cuts, we included data up to 0.1 dcx to the 
right of the diagnostic line to account for possible uncer- 
tainties. 

Roughly 25% of our spectra were rejected because of 
poorly-fit H7 emission, line ratio cuts, or inadequate line 
fits. Spectra rejected due to line ratio cuts are shown as 
black crosses atop the galaxy images in Figure |4j Our 
final sample of H II regions in early-stage interacting 
galaxies consists of 332 spectra. These are accompanied 
by 281 H II region spectra from the control sample. Four 
nuclei in our sample show either Seyfert (NGC 2207, IC 
2163, NGC 7469) or LINER (NGC 3994) spectra based 
on our data. However, these nuclear ionizing sources 
have minimal impact on the extra-nuclear H II regions 
that are used to compute abundances; phenomena like 
ionization cones would in any case be detected and re- 
jected using our line ratio cuts. 

In this paper, we quote abundances calculated from 
the [N II1A6583/fO II1AA3727 , 3729 abundance diagnos- 
tic ( jKewley fc Dopita||2lj02"} hereafter KD02). Of the 
strong- line diagnostics, this one produces the least root- 
mean-square (RMS) dispersion about the mean due to 
its relative insensitivity to ionization parameter (KD02) . 
To study possible systematic errors from the use of this 
diagnostic (due to, e.g., improper extinction corrections 
or improper relative flux calibration of the blue and red 
LRIS arms), we compu ted abundances using the Kobul-| 



nicky fe Ke wley (2004, hereafter KK04) R 23 and Pettini 
fc Pagel||20(jl^ eTeairer PP04) ([O III]/H/3)/([N Il l/Ha) 



diagnostics (follow ing the guidelines laid out by |Kew 
ley fe Ellison|2008 ) and then converted these dia gnostics 
into the KL)(J2 c alibration using the formulas of jKewleyl 
& Ellison (2008). The results are shown for all of the 



ti 11 regions m this paper in Figure [3j and for individual 
systems in Figure [4] 

It is clear from examining Figure [3] that, overall, 
there are no significant systematic discrepancies intro- 
duced by our data reduction and analysis that make our 
[N II]/[0 II] line ratios suspect. For instance, both KK04 
and PP04 are less sensitive to extinction than KD02, 
and inaccurate extinction corrections should show up as 
deviations from equality in these figures. The interact- 
ing galaxies sample does in fact show a 0.05 dex offset 
from equality in the KD02 vs. PP04 diagram. How- 



ever, the PP04 diagnostic does not correct for ionization 
parameter variations, which introduces a systematic un- 
certainty in this diagnostic; the KD02 diagnostic is insen- 
sitive to ionization parameter. Furthermore, the offset is 
well within the expected uncertainties of the diagnostic 
and conversion formulas, and as a simple offset, will have 
no impact on relative abundances like gradients. In fact, 
the RMS dispersion about the line of equality in the in- 
teracting galaxies sample (0.05 dex) is significantly lower 
than that in the control sample for the KD02 vs. KK04 
diagram, suggesting that our data is highly uniform. In 
the KD02 vs. PP04 diagram, the RMS dispersions are 
the same for the two samples, but for the interacting 
sample this is in large part due to the systematic offset 
to one side of equality. 

Errors in the emission-line fluxes were computed by 
assigning the error in the peak flux to equal the RMS 
residual within 2a of each line center, after subtracting 
the continuum and line fits from the data. Errors were 
propagated primarily using analytic expressions, but for 
abundance and gradient errors we employed Monte Carlo 
methods. 

In the following section, we use the calculated abun- 
dances and radii for the H II regions in each galaxy to 
compute radial abundance gradients, and compare the 
gradients of the control and interacting samples. 

4. OXYGEN GRADIENTS 

With abundances and deprojected radii in hand, we 
computed radial oxygen gradients for each galaxy in our 
sample. In two cases (IC 5283 and UGC 12915), the ra- 
dial coverage and number of H II regions available were 
each too small to yield a reliable estimate of the gradi- 
ent. In the control sample, we limited the gradient fits 
to within I.5R25, to match the radii fit in the interacting 
galaxy sample. We wished to avoid biasing the gradi- 
ent fit toward particular locations in the galaxy (e.g., 
to high- luminosity H II regions), so we performed un- 
weighted least-squares fits. Our fits were computed with 
radius as the independent variable, and errors in slope 
and intercept were estimated using Monte Carlo meth- 
ods. The resulting gradients and intercepts are given in 
Table [p 

In Figure [4l we present for each system the results of 
our emission-line analysis: emission line diagnostic dia- 
grams, comparison of different metallicity diagnostics, a 
two-dimensional projected oxygen abundance map, and 
a one-dimensional deprojected radial oxygen abundance 
profile. 

The measured gradients in these interacting systems 
are on average significantly flatter than those in the con- 
trol sample (Figure pi). The control sample exhibits 
median slopes of -0.0Tl±0.009 dex/kpc and -0.57±0.05 
dex/R 2 5, while the interacting systems have median 
slopes of -0.017±0.002 dex/kpc and -0.23±0.03 dex/R 25 
(where the errors given are the standard error). In short, 
the interacting systems have gradients that are, on aver- 
age, less than half as stee p as isolated galaxies. 

we showed that gradients in 
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eight of the systems presented here were shallower than 
gradients in M83, M101, and the Milky Way. The present 
data quantify this difference at a high confidence level 
using large control and interacting samples. Taken by 
themselves, these differences in gradients are a strong 
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confirmation of the model of merger-induced gas inflow 
leading to metal mixing. 

Surprisingly, the dispersions in the interacting sample 
gradient distributions are also less than the dispersions 
in the control sample: the control sample shows standard 
deviations of 0.03 dex/kpc and 0.17 dex/R 2 5, compared 
to 0.01 dex/kpc and 0.11 dex/R,25. The explanation for 
the latter effect is unclear, but could result from the ori- 
gin of the H II region data in heterogeneous and homo- 
geneous datasets, respectively. 

No significant correlations are seen between gradients 
and galaxy near infrared and total infrared luminosities, 
either in the interacting sample or the control sample 
(Figure [7]) . Previous work has shown that correlations 
may exis t between optical lu minosities and gradients in 
dex/kpc ( Zaritsky et al.|1994[) ; however, such a trend ap- 
parently does not persist m interacting systems (and is 
not obvious in our control sample). One might expect a 
correlation of infrared luminosity with gradient in inter- 
acting systems, if Lir is assumed to trace star formation 
rate and star formation is driven by gas inflow. Again, no 
such trend is evident, though we only have upper limits 
in Lyn for much of our sample. 

No obvious trends are seen when comparing gradients 
to system properties that may parameterize the stage or 
strength of the interaction: projected nuclear separation, 
near-infrared luminosity ratio, or total system infrared 
luminosity (Figure [8]). Projected nuclear separation is of 
course only approximate, since we cannot deproject the 
two galaxies without detailed modeling, and there is a 
degeneracy in the age of the interaction between galax- 
ies moving apart after first pericenter and returning to 
second pericenter. Furthermore, mass ratios in our sam- 
ple (as traced by -ftf-band luminosity ratio) are all close 
to unity, and one might expect little variation among the 
systems in our sample based on this parameter alone. Fi- 
nally, confusing the matter, simulations also show that 
although both galaxies in a strongly interacting pair ex- 
perience a flattening of the gradient very quickly after 
first pericenter, the relative degree of flattening of the 
two disks varies dependi ng on the geometry and the dis- 
tance of first pericenter ( Rupke et al.||2010 ). This effect 
would introduce scatter mto any underlying cor relat ion 
between gradients and system properties. See S5.2 for 
further discussion. 

Some systematic uncertainty arises in our abun- 
dance gradients due to our choice of abundance cali- 
bration/diagnostic. Some calibrations (including both 
strong- and weak-line ones) may yield shallower slopes 
than we calculate (e.g. JKewley fc Ellison 2008 Bresolin 



et al.|200"9a ). However, kewlcy et al. (2010]) showed that 



an R 2 3- based diagnostic yields the same results as the 
KD02 diagnostic we employ in this paper. More impor 
tantly, we rely on the fact that within a particular abun 
dance calibration/diagnostic, relative abun dances are re 



liably com puted and internally consistent ( Kewley & El- 
lison||2008|). When using detailed chemodynamical mod- 
els to constrain parameters of individual systems, mea- 
suring the absolute value of the gradient precisely may 
be important. Within the context of the current study, 
we are concerned only with gradient changes relative to 
initial conditions. We also emphasize that, were the ini- 
tial galaxy gradients too shallow (or completely flat), 
the metal dilution clearly observed in interacting systems 



would be impossible, and would thus be inconsistent with 
empirical evidence. 

Nevertheless, as a check on the reliability of our 
[N II] /[O II] gradients, Figure [1] shows, for each system, 
how the [N II] /[O II] abundances compare with t hose 
computed in the KK04 and PP04 diagnostics (see ^3.2 
for more details). As is seen in the full sample (Fig. [3]), 
individual systems are largely free of systematic discrep- 
ancies. The two exceptions are NGC 2207/IC 2163 and 
Arp 248, which show discrepancies between the KD02- 
and PP04-computed metallicities. These anomalies may 
result from variations in ionization param eter that re- 
main uncorrected in the PP04 diagnostic 3.2 ) or red- 
dening uncertainties. Further analysis shows that, if we 
were to use the PP04 abundances (after being converted 
into the KD02 system) , we would measure shallower gra- 
dients in Arp 248 NED01 and Arp 248 NED02 and a 
steeper gradient in NGC 2207 (but not IC 2163), in each 
case by a factor of ~2. 

With the common exception of the galaxy nuclei, the 
vast majority of data in our sample are pure H II regions. 
However, UGC 12914/12915 has a significant number of 
regions indicative of non-photoionized emission. In fact, 
the north half of UGC 12914 has a significant number 
of regions with multiple velocity components and shock- 
like line ratios ( |Allen et al |2008| |Rich et al.||2010| . This 
region overlaps a bri dge ot H I gas between the galaxies 
( |Condon et at|p93l ). 



5. COMPARISON TO SIMULATIONS 

Numerical simulations are a powerful tool for investi- 
gating the evolution of galaxy mergers. However, most 
simulations of mergers to date have paid little atten- 
tion to gas-phase m etallicity and its evolution. Recently, 
Rupke et al. (2010) discussed metallicity evolution prior 
to second pericenter in smoothed particle hydrodynamic 
(SPH) / N-body simulations of equal-mass mergers with 
gas mass fractions of 10% (to match nearby galaxies). 
Though they d id not include enric hment from ongoing 
star formation, |Rupke et al.| ( |2010|) reached good agree- 
ment with the magnitude ot deviations of major mergers 
from the mass-metallicity relation. They also showed 
that radial abundance gradients should flatten rapidly 
after first pericenter, and predict the relationship be- 
tween gradient flattening and nuclear abun dance and 
the d e tailed shape of the gradient with time. Montuori 
(2010[) presented simulat ions including star fb rma- 
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tion that largely agreed with |Rupke et aL (2010), but 
did not address the evolution ot gradien ts, lorrey et al. 
(2010, in prep.) find similar results to Montuori et al. 
(20101, but show that preferential consumption ot high- 
met allicity gas by star formation in galaxy nuclei plays 
a role alongside metal dilution due to inflow. 

In the following sections, we compare our data on 
metallicity di stributions in intera cting galaxies to the 
simulations of |Rupke et al | ( |2010j ). These SPH/N-b ody 
simulations are described in detail in Barnes (2004 1. In 
brief, 24576 gas particles were distributed in an expo- 
nential disk of scale height equal to 6% of the disk scale 
length. The stars occupied a bulge and identical disk, 
with a bulge-to-disk ratio of 1/3. Eight close-passage 
mergers were simulated, with all combinations of direct 
and retrograde considered. The metallicity per gas par- 
ticle was unchanged during the simulation, and two dif- 
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ferent initial metallicity gradients were considered (-0.2 
dex/Rdisk and -0.4 dex/Rd). In this work we modified 
the initial gradients slightly to match our control sam- 
ple: -0.15 dex/Rd- This number results from the typical 
gradient in dex/R 2 5 for the control sample (Q, when 
c onve rted to dex/Rd using a typical ratio for~R25/-Rd 

(S 

5.1. Gradient Distribution 

Our galaxy sample was chosen to represent major 
mergers somewhere between first and second pericenter. 
The distribution of observed gradients is shown in Fig- 
ure [5j For comparison, we construct a predicted gradient 
distribution from our set of simulated major mergers, as- 
suming that our sample is randomly distributed between 
first and second pericenter. We also assumed that the 
simulated galaxies start with a single gradient equal to 
the mean value in the control sample. The observed and 
simulated distribution show remarkable agreement (Fig- 
ure [9}. 

5.2. Gradients vs. Merger Stage 

We show in Figure [8] that, for the present sample, there 
is no strong dependence of gradient on merger stage, as 
parameterized by projected nuclear separation. This is 
consistent with the predictions of our simulations. Fig- 
ure 10 reveals that, between first and second pericenter, 



We st a rt wit h the R,23-computed abundances from|Salzer| 



et al. (2005); the diagnostic they use is a functional fit 
of K.23 to abundances from modeling of SPSS sp ectra 
(|Charlot fc Longhetti||2001| |Tremonti et al.||2004). We 



then apply the formula from |Kewley &; Ellison |2008 ) . 
Because two transformations are involved, we consider 
these abundances to be only approximate, and are useful 
primarily for establishing the slope of the L — Z relation. 

We do not have nuclear abundances measured from 
long-slit spectra for all of the galaxies in the current sam- 
ple. We thus estimate abundances by extrapolating our 
gradient fits, and compute the abundances at a discrete 
radius of O.IR25. This ro ughly matches the extent of the 
long-slit apertures of the Salzer et al. ( 2005 ) data, which 
have (z) = 0.063. 

Using this method, neither the control nor interact- 
ing samples fall on the i^-band L — Z relation. This 
is pre sumably due to one of two effects: (1) the Salzer 



et al. (2005) abundances are luminosity-weighted mea- 
surements from nuclear spectra, while ours are measured 
at a fiducial radius using fits to H II region data; and/or 
(2) the uncertainties introduced by the two trans forma- 
tions involve d in computing the abundances of thejSalzer 



there is no predicted dependence of the gradients on ac- 
tual nuclear separation, for the ensemble of simulations 
as a whole. Within particular mergers, a more recogniz- 
able dependence of gradient on nuclear separation may 
emerge, but we must consider our present dataset as ran- 
dom snapshots of an ensemble of different mergers. We 
do not at present have an accurate understanding of the 
initial conditions of most or all of the systems in our 
sample. 



et al.| (I2005I) sample. The latter effect is not unusual 
when converting among different metallicity calibrations 
(Kewley & Ellison [2008 ) . To account for this effect, we 
apply a constant offset of —0.3 dex to both the control 
sample and interacting sample metallicities. This offset 
minimizes the RMS deviation of the control sample from 
the L — Z relation. 

The exact vertical normalization of the data is irrele- 
vant, however; the important quantities are the relative 
nuclear abundances of the control and interacting sam- 
ples at a given luminosity, which are secure since we use 
the same abundance diagnostic and compute the nuclear 
abundances in the same way. The Salzer et al. (2005) 
data simply provide a more accurate slope tor the L — Z 



5.3. Gradients and the Lummosity-Metallicity Relation relation than the control sample by itse lf. 



The motivation for the present work was the observa- 
tion that interacting galaxies fall below the luminosity- 
metallicity (L—Z) and mass-metallicity (M—Z) relations 
of star-forming galaxies. Do the current data present a 
consistent picture in this regard, such that the nuclear 
metallicities of interacting systems with shallow gradi- 
ents also fall below these relationships? For a subset of 
th e current sample, w e used the S-band L — Z relation 



The result of this exercise is displayed in Figure 11 



Kewley et al. (2010) and found no correlation between 



gradients and L — Z offsets, but were limited by small 
sample size and the fact that -B-band luminosity is sen- 
sitive to ongoing star formation and extinction effects. 

In this paper, we compare the control and interact- 
ing samples to the luminosit y-metallicity relatio n as mea- 



sured in the near-infrared (Salzer et al 
simulations to interpret the result. 



2005 ), and use 
Multicolor broad- 



band photometry is not available for our entire sample, 
so we are unable to measure accurate masses and apply 
the M — Z relation. Near-infrared photometry, however, 
is less sensitive to extinction or contamination from star 
formation than the optical; it is also a reasonable tra cer 
of stellar mass within a factor of ~2 ( Marasto"n"||2005 1 . 

The Salzer et al. (20051 abundances are available in 
several diagnostics. None of these are the [N II] /[O II] 
diagnostic that we employ, s o we resort to the abun- 
dance conversion formulas of Kewley & Ellison ( 2008 ) . 



This figure makes evident that the early-stage interacting 
galaxies, which have shallower gradients than the control 
sample, also fall below the L — Z relation defined by the 
control sample and emission line galaxy sample. In fact, 
the interacting galaxies form a remarkably coherent L—Z 
relationship that is simply offset from that of isolated 
star-formers by —0.2 dex in abundance. 

To tie together abundance offsets and gradients, we 
compa re directly to the simulations of |Rupke et al.| 
(2010). These simulations predicted that the gradient 
slope would change quickly after first pericenter, while 
the nuclear metallicity would change more slowly be- 
tween first and second pericenter. The simulations yield 
gradients in units of exponential disk scale length of the 
gas disk (dex/R d ), while we express observed gradients 
in physical units (dex/kpc) and optical (stellar) isopho- 
tal radius (dex/R2s). For our control sample, we have 
stellar disk scale lengths as well as optical radii, and the 
two are related according to R25 = (3.7 ± 1.2)Rd, where 
the uncertainty given is the standard deviation. Using 
this scaling, we compute the gradients in dex/R^ for our 
interacting galaxy sample. 

For the current data set, we measured abundances at a 
fiducial radius of O.IR25. To match this procedure for the 
simulations, we computed average abundances within a 
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narrow, O.lRrf-wide radial bin centered on 0. 35Rrf. (This 
differs from the method in Rupke et al.|2010| in which we 
averaged over the entire central disk within 0.5R<2.) For 
the data, we assume that the interacting galaxies started 
on the L — Z relation of the control sample and have had 
their central metallicity lowered during the merger. 

It is clear that, given the uncertainties inherent in this 
procedure (e.g., the galaxies didn't necessarily start on 
the mean L — Z relation) and the simplicity of the sim- 
ulations (which do not include ongoing star formation), 
there is remarkab le a gree men t between the data and sim- 
ulations (Figures fl2| and [T3l>. A total of 70 - 80% of the 
interacting galaxies fall directly in the phase space de- 
lineated by the simulated evolutionary tracks. Further- 
more, the region of overlap with the simulations is near 
the time of first turnaround (between first and second 
pericenter) for many of the simulated pairs, which is con- 
sistent with the galaxies in our sample being somewhere 
between first and second pericenter on average. (In a 
random distribution of post-pericenter pairs selected by 
nuclear separation, half should be pre-turnaround, and 
half post-turnaround.) 

The data are thus quantitati vely consistent wit h the 
prediction of the simulations of |Rupke et al.| ( 2010 ) that 
the gradient should flatten quickly after first pericenter, 
and that the nuclear abundance should decrease at a 
slower rate between first and second pericenter. 

The scatter of plausible initial conditions, as repre- 
sented by the control sample, do caution against drawing 
detailed conclusions from this comparison of data and 
simulations. In particular, a variety of initial conditions 
can explain the current position of the interacting galax- 
ies in the metallicity-gradient phase space. The metallic- 
ity scatter in the control sample is partly a measurement 
issue; e.g., the scatter in the simulated control L — Z 
relation is larger than in the Salzer et al. (20051 sam- 
ple, and the control sample H II regions show a larger 
RMS dispersion in abundance than the interacting sam- 
ple (Fig. [6]). To deal with this metallicity scatter, what- 
ever its origin, howeve r, w e simply shift the simulations 
left or right in Figure [12] to start at the control sample 
initial conditions. 

A shift in L — Z offset could in principle resolve any 
discrepancies between data and simulations for individ- 
ual galaxies. For example, NGC 2207 is clearly not at a 
late merger stage, as suggested by its position with re- 
spect to the merger models (it is probably just after first 
passage; Elmegreen et al. 1995[ ). However, it could ini- 
tially have had lower metallicity than predicted by the 
mean L — Z relation and/or a steeper initial gradient, 
leading to a more favorable comparison with models in 
Figure [T2] 

If a galaxy were to start with a larger gradient than is 
assumed in the models, then t he changes in gradien t and 
central metallicity are larger (Rupke et al. 2010). For 
instance, if we assume the initial conditions represented 
by NGC 300 or NGC 2997, the simulations still overlap 
the data, due to the more extreme changes that result 
from a stee per gradient and the same amount of infall 
(Figure 14 1. Closer inspection rules out this scenario for 
most systems, however, since (a) the overlap of the data 
and the bulk of the simulations is not as good and (b) 
these initial conditions would require almost all of the ob- 
served interacting systems to lie near second pericenter. 



T his is unlikely, due to how we have selected our sample 
(j ]2.1|). The simulation-to-data comparison shown in Fig- 
ure 1 12 [ with half the initial gradient of that in NGC 300 
or NGC 2997, is therefore a more likely scenario for most 
systems. 

5.4. Gradient Profiles 

Following standard practice, we have characterized 
radial metallicity profiles using straight line fits [i.e., 
12+log(0/H) = a + bR]. However, t he predicted metal- 
licity profiles from the simulations of Rupke et al. (20101 
do show distinctive shapes that are not purely straight 
lines at all times. 

Individual galaxies in our sample do not have enough 
data points to accurately characterize the profile shapes 
beyond straight lines. However, we can compute the av- 
erage deviation from straightness as a function of radius 
by subtracting the gradients and binning the data in ra- 
dius. We do this for both the control and interacting 
samples, and show the results in Figure |15| These two 
samples show metallicity profile residuals that are signifi- 
cantly different. The control sample is largely consistent 
with a straight line. However, the interacting sample 
shows a clear "concave up" shape, in the sense that pos- 
itive deviations are seen at small and large radius, with 
negative deviations at middle radii. 

How does this compare with simulations? To bin the 
simulation data, we consider the ensemble of profiles at 
first and second pericenter and first turnaround. The 
result, with the interacting galaxies data overplotted, is 
shown in Figure 15 It is apparent that for the sim- 



ulations, there is a small upturn at low radius at all 
times, consistent with the data. Though there is not 
exact quantitative agreement in any of the three stages, 
the simulations at first turnaround show the same "con- 
cave up" shape seen in the data. By contrast, the sim- 
ulations at 1st and 2nd pericenter show a flatter profile 
outside the nuclear regions. As discussed above, this is 
consistent with the idea that, on average, our sample is 
near first turnaround, and otherwise evenly distributed 
between first and second pericenter. 

Given the complexity and variety of the actual merging 
systems in this study, it is impressive that there remains 
such good agreement at this level of comparison between 
data and simulations, and speaks to the robustness of 
the dilution model and its predictions. It also suggests 
that metals produced in star formation initiated by the 
merger have only a small effect on the gas-phase metallic- 
ity distributions prior to second pericenter. This could 
occur either because most star formation occurs after 
second pericenter or because of inefficient mixing of the 
metal-enriched hot phase with the cooler ISM. In fact, 
only a limited fraction of the gas in interacting galaxies 
can be consumed in star formation prior to merging, in 
order to leave a sign ificant gas reservoir at late stages 
(Sanders et al.|1991|. 



6. SUMMARY AND PROSPECTS 

We have presented here data on the gas-phase oxygen 
abundance gradients in 9 pairs/groups (18 galaxies) in 
the early stages of an interaction (after first pericenter, 
prior to second). We computed radial gradients based 
on H II region spectroscopy and the most precise strong- 
line abundance diagnostic available (KD02). We find a 
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clear break between the gradient distributions of isolated 
and interacting systems: interacting systems consistently 
have gradients about half that in normal, isolated spirals. 

We comp are the data to the simulations of [Rupkc 
et al. (2010), and find remarkable statistical agreement 
at the ensemble level. The simulations accurately repro- 
duce (1) the gradient distribution, (2) the lack of gradi- 
ent dependence on projected nuclear separation, (3) the 
relationship between gradients and deviations from the 
luminosity-metallicity relation, and (4) metallicity profile 
deviations from a straight line. 

These data are clear evidence that models of metal re- 



distribution due to interaction-induced gas inflow |Kew- 



ley et al.|2006a||Rupke et al.|2010| [Montuori et al.|2010| ) 

are able to account for the low nuclear abundances ob- 



served in interacting galaxies ( Lee et aI7||2004[ |Kewley| 
et al. |2006a[ |Rupke et al. |2008[ |Kllison et al. I^WTa) |Michel- 
Dansac et al.|2008||Peeples et al.|2009||Ekta <E"C hcngalur 
2010||Sol Alonso et al.||2010p . Now that this mechanism 
is firmly established, future work should be devoted to 
understanding how this affects the global evolution of the 
mass-metallicity relation over cosmic times and the for- 
mation of elliptical galaxies. In particular, if mergers or 
post-mergers dominate samples of galaxies used to deter- 
mine the M — Z relation at high redshift, this bias must 
be taken into account. Furthermore, all massive galaxies 
have experienced mergers during their lifetime. The re- 
distribution of metals due to gas inflow will then affect 
the stellar populations formed at later times at all radii. 

Though we have already found remarkable agreement 
between the data and simulations at a statistical level, 
this comparison could be improved upon using simu- 
lations that incorporate realistic models of star forma- 
tion and chemical enrichment (including a sophisticated 
model of starburst-driven galactic winds). Secondly, we 
have compared to simulations only after binning the data 
azimuthally; more detailed comparison of theory to ac- 
tual projected metallicity distributions would provide a 
stronger understanding of how and where the dilution 
occurs. Furthermore, we have only looked for agreement 
at a statistical level. Though the current data already 
provide broad constraints on the simulations, more de- 
tailed constraints will require careful modeling of indi- 
vidual systems. 

The limited merger ages traced by this early-stage sam- 
ple constrain our ability to trace metal redistribution 
over the entire course of a merger. In future work (Rupkc 
et al. 2010, in prep.), we will present data on later stage 
mergers to accomplish this. 
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Figure 1. Properties of the control and interacting galaxies samples: (a) K-band absolute magnitude; (b) 8— 1000/im infrared luminosity; 
(c) iC-band luminosity ratio; and (d) nuclear separation. The control and interacting samples are shown as red and blue hatched histograms. 
In (b), the open histogram is the total system luminosity for the interacting sample. In (c), the convention is such that the luminosity ratio 
is always >1. The control and interacting samples are very similar in Mjf, but the interacting galaxies have higher infrared luminosities. 
The interacting galaxies have luminosities near L* , mass ratios in the range 1:1—1:2.5, and (projected) nuclear separations in the range 
15 - 30 kpc. 
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Figure 2. Example fits to a high-flux and low-flux H II region spectrum from our sample (the latter has 1% of the Ha flux of the former). These demonstrate the high 
fidelity of our fits over a large range of signal-to-noise. The top panels show close-ups of various emission line fits, with the residuals underneath each fit. The bottom 
panels show the fits to the stellar continuum, with the emission-line fits subtracted. Black is our data, and red is the fit. The spectra are shown in the observed frame. 
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Figure 2. Continued. 
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Figure 3. Comparisons between different strong-fine methods of cafcufating abundances, within the control sample (red open circles) 
and the interacting sample (blue fitfed circies). The diagonai bfack line shows equaiity, and the biue and red dashed fines show the RMS 
dispersion o f each sampfe (0.05 — 0.1 de x) ab out equality. We compare our chosen diagnostic (KD02) with KK04 and PP04, converted into 
KD02 using [Kewley fc Ellison] J20 08 > (j ]3.2| l. A smalf shift (0.05 dex) is seen toward higher metaificity when comparing KD02 and PP04, 
but this is well withm the uncertainties of the calibrations and cross-calibrations. The good agreement between these methods points to 
the absence of strong systematic errors in our measurement of line fluxes. 
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Figure 4. Results of emission-line analysis of each system, including a map of the oxygen abundance from the [N II]/[0 II] method, ovcrplotted on an Ha image (or 
F435W HST image in the cases of Arp 256 and Arp 298); abundance gradients for each galaxy, versus radius in kpc and radius normal ized to R25; comparisons between 
the KD02 abundance diagnostic and the KK04 and PP04 diagnostics; and the emission line ratio diagnostics from |Kewley et a!"] J2006b) . In each plot, the data points arc 
colored according to the KD02 abundance (sec color bar). The images show the near-infrared nuclear position as a white open circle, and a horizontal bar in the upper-left 
corner represents 5 kpc. The gradient panels show our formal fit and t he average RMS (0.05 de x) around this. For the abundance comparisons, the KK04 and PP04 
abundances a re converted to the KD02 diagnostic using the conversions of Kewley &: Ellison 1 2008 I. We overplot the line of equality and the RMS scatter in the |Kewley &| 
|Ellison|(2008) SDSS sample after this conversion (~0.05 dex). In the emission line diagnostics, points used in the abundance calculations are solid symbols; poinis rejected 
based on the [N II]/Hct and [O I]/Hct diagrams are open symbols. These rejected points are plotted as black crosses on the images. 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 4. Continued. 
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Figure 6. Distribution of RMS residuals from a straight-line gradient, for the interacting and control samples. The interacting sample 
shows smaller, rather than larger, RMS residuals. This implies that the metallicity scatter in galaxies does not increase due to a merger; 
the observed decrease in the scatter is likely an artifact of the relative homogeneity of our data compared to the control sample. 
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Figure 7. Gradients vs. galaxy properties: iC-band and total infrared luminosities. Open circles are the control sample, and closed circles 
are the interacting sample. Arrows represent systems with unresolved infrared luminosity; in these cases the total system luminosity is an 
upper limit to the luminosity of each galaxy. No significant trends are apparent within either subsample. 
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Figure 8. Gradients vs. system properties: nuclear separation, mass ratio (as traced by X-band luminosity), and total infrared luminosity. 
Pairs are connected by straight lines, and the most massive galaxy in each pair is circled. No significant trends are apparent. 
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Figur e 9. Distribution of gradients in the interacting sample compared to the distribution predicted from the simulations of Rupkc ct al. 
( |2010| l. The blue, dashed vertical line shows the average gradient in the control sample, and the red, solid vertical line shows the initial 
gradient for each simulation. The simulated distribution is computed assuming that the galaxies lie between first and second pericenter. 
Very good agreement is evident. 
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Figure 10. Simulated dependence of gradient on (actual) nuclear separation between first and second pericenter jRupke et al.||2010| . 
Each point represents a time step for one of two galaxies in one of eight simulations. Points are connected for each galaxy, witn a different 
color assigned to each simulation. The horizontal solid line is the initial gradient for each simulation. At the statistical level, no clear 
dependence of gradient on nuclear separation is expected, consistent wit h what is observed (Figure jg}. It is also evident that in most 
systems, the gradient flattens quickly after first pericenter, as discussed in |Rupke et al.| | |2010| l. 
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Figure 11. Near-infrared luminosity-metallicity relation. The small black points are emission-line galaxies from [Salzer et ah] ( |2005| ; the 
open red circles are the control sample; and the filled blue circles are the interacting galaxies. Nuclear abundances tor our sample arc 
computed at R = 0.1R25- The solid line is the bisector of the Y on X and X on Y unweighted least-squares fits to the |Salzer et al.] ( |2005} 
sample, with RMS dispersion illustrated by the dashed lines. The control and interacting samples are shifted down by 0.3 dex so that the 
control sample lies atop the L — Z relation. The interacting sample falls well beneath the L — Z relation delineated by the control and 
Salzer ct al. (2005) samples. 
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Figure 12. Estimated offset from the L — Z relation vs. metallicity gradient. Open circles are the control sample, and filled cir cles are 
the interacting sample. Dashed blue lines illustrate the time evolution in this phase space of the eight major merger models from [Rupke] 
|et al.| | |2010| . The models move from the center to the upper left. The triangles, squares, and star indicate various important points in the 
evolution or the simulated merger. The input gradient for the simulations is chosen to match the average gradient per stellar disk scale 
length from the control sample. Given both the measurement and systematic uncertainties in this procedure, there is excellent agreement 
between the data and models. The data are consistent with the typical system lying near first turnaround. 
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Figure 13. Estimated offset from the L — Z relation vs. metallicity gradient. See Figure [12] for more details, fndividual galaxies are 
indicated. 
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Figure 14. Estimated offset from the L — Z relation vs. metallicity gradient. See Figure [T2] for more details. In the simulations displayed 
here, we used for initial conditions the steepest gradients observ ed i n the control sample, and the corresponding L — Z relation offsets. 
The overlap between data and simulations is poorer than in Fig. |12| and this scenario would require most of the pai rs to lie near second 
pericenter rather than first turnaround. The data thus favor the shallower progenitor gradients illustrated in Fig. |12| 
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Figure 15. (left) Metallicity residual from a straight line fit vs. radius for all H II regions in the control (open red circles) and the 
interacting (filled blue circles) samples. Binned points are averages over bins of size O.25R25, with the standard error shown as vertical 
error bars. The control sample is consistent with a mostly flat profile, while the interacting sample shows a "concave up" shape, (right) 
The same, but for the simulations (black points and red, open circles), with the binned data for the interacting sample overplotted. The 
best agreement in shape is found between the interacting sample and the simulations at first turnaround, which both have a "concave up" 
profile residual. Exact quantitative agreement is lacking, but this comparison is consistent with the idea that our sample lies on average 
near first turnaround. 
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Table 1 

Sample 



Galaxy 


System 


cz 


Nuc. Sep. 


M K 


kWLrn/Lp,) 

O \ 1-1- W \ZJ / 


R25 


i 




Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) ES 


(10) 


UGC 12914 


VV254 


4371 




-24.23 


10.36 


14.69 


54 


160 


2,3,4 


UGC 12915 


VV254 


4336 


17.40 


-23.70 


10.76 


11.11 


73 


135 


2,3,4 


Arp zoo iNH/JJUz 


Arp 256 


8193 


' ' ' 


-23.11 




16.58 


0—45 


t u 


1,3,4 


Arp ZOO iNH/JJUl 


Arp 256 


8125 


28.20 


-23.80 


11.35 


13.66 


61 


80 


1,3,4 


UCjU olz 


HCG 2 


4364 




-21.56 


9.93 


11.67 


79 


8 


2 


Mrk ooz 


HCG 2 


4349 


21.10 


-22.33 


10.38 








2 


T Qi a a 


HCG 2 


4271 


65.40 


-20.52 












UGC 813 


WBL 036 


5344 




-22.89 


10.68 b 


11.16 


72 


110 


2 


UGC 816 


WBL 036 


5188 


16.50 


-23.01 


10.68 b 


13.33 


62 


190 


2 


CGCG 551-011 


WBL 036 


5373 


42.60 


-23.27 










2 


1NCC 2207 


NGC 2207/IC 2163 2741 




-24.54 


11.00 


22.11 


35 


140 


3,4 


IC 2163 


NGC 2207/IC 2163 2765 


15.80 


-24.18 




9.80 


35 


128 


3,4 


Arp 248 NED01 


Arp 248 


5167 




-22.39 


10.66 b 


13.90 


66 


43 


1 


Arp 248 NED02 


Arp 248 


5167 


54.60 


-23.32 


10.66 b 


13.27 


0-45 


109 


1 


Arp 248 NED03 a 


Arp 248 


5276 


93.10 


-20.82 










1 


NGC 3994 


Arp 313 


3086 




-23.56 


10.55 


6.18 


58 


10 


1.2,3,4 


NGC 3995 


Arp 313 


3254 


24.70 


-22.64 


10.42 


19.40 


69 


45 


1,2,3,4 


NGC 3991 a 


Arp 313 


3192 


50.10 


-22.51 


10.26 








4 


NGC 7469 


Arp 298 


4892 




-24.23 


11.46 


12.22 


30 


127 


1,3,4 


IC 5283 


Arp 298 


4804 


23.30 


-23.36 




9.90 


60 


105 


1,3,4 


UGC 12545 


WBL 713 


5754 




-21.41 


10.28 b 


11.81 


63 


85 


2 


UGC 12546 


WBL 713 


6055 


21.90 


-22.69 


10.28 b 


10.64 


63 


20 


2 



References. — 1 = |Arp| ( fl966t >; 2 = |Barton et a.l.| ( |2000| ); 3 = |Sanders et a.l.| ( |2003] ) ; 4 = |Surace et al.| l |2004[ | 
Not e. — C ol.(l): Galaxy name. Col. (2): System label from optical identificati on. In order of preference : Arp 



1066 



VV = |Vorontsov-Velyaminov| (1977) ; HCG = Hickson Compact Group =[ffickson (1982]; WBL = |White et 



19991 catalog of poor clusters. Col. (3): Rcdshift from NED, in km s . Col. (4): Projected nuclear separation from 
first galaxy listed for system, in kpc, using. Col. (5): A*-band ma gnitude down to 20 mag/a rcsec 2 isophote, from 2MASS. 
Col. (6 ): 8— 1000/^m luminosity, determine d using the formula in |Sanders &i Mirabel| | |1996[l and fluxes from |Surace et ak] 
( |2004| > , the IRAS Faint Source Catalog, or |Johnson et ah |||2007^ in the c ase ot riUU 2. Uol. (V ) : B-band optical radius trom 
25 mag/arcsec 2 isophote, in kpc, from HyperLeda |PatureTet al.|2003| ). Col. (8): Inclination in degrees, from HyperLeda 
except for Arp 248 NED02 and Arp 256 NED02. For these cases, we assume 30°; the actual inclination is in the range 
— 45°. Col. (9): Position angle of the galaxy line of no des, in degrees east o f north; from HyperLeda except for Ar p 25 6 
and WBL 036 (used H I data from |Chen et al.|2002| and |Condon et al.|2002| ). Col.(10): Sample selection reference (§231. 
a No LRIS data is available for these galaxies. They are listed here as prominent but more distant members of the relevant 
systems (except for Mrk 552, which is the closest galaxy to UGC 312 in projection). 

b The infrared luminosity for this pair is unresolved; the total system luminosity is computed using the redshift of each 
galaxy. 
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Table 2 

Control Sample 



Galaxy 


Type 


Distance 


Dist. Ref. 


M K 


log(£iR./L©) 


R.25 




R d Ref. 


i 


PA 


i/PA Ref. 


H II Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 300 


SAd 


2.08 


a 


-19.43 


8.43 


5.90 


2.26 


B 


40 


114 


aa 


1 


NGC 628 (M74) 


SAc 


8.59 


b 


-22.48 


9.82 


13.08 


3.57 


A 


7 


25 


cc 


6,7 


NGC 925 


SABd 


7.38 


c 


-20.75 


9.21 


12.04 


2.60 


D 


61 


107 


aa 


7 


NGC 1232 


SABc 


14.50 


d 


-23.08 


9.91 


14.59 


4.51 


C 


33 


270 


<\'<\. CO 


3,7 


NGC 2403 


SABcd 


3.16 


a 


-21.05 


9.19 


9.17 


1.72 


C 


60 


127 


aa 


1,7 


NGC 2805 


SABd 


28.00 


c 


-21.44 


9.68 


14.42 


8.96 


E 


36 


290 


aa,bb 


7 


NGC 2903 


SABbc 


8.55 


e 


-23.56 


10.22 


14.95 


3.09 


E 


56 


23 


aa 


3,7 


NGC 2997 


SABc 


7.08 


c 


-22.61 


9.82 


11.55 


4.84 


F 


41 


110 


aa 


3 


NGC 3184 


SABcd 


8.70 


c 


-22.09 


9.48 


9.60 


2.59 


G 


24 


135 


aa 


7 


NGC 5236 (M83) 


SABc 


4.92 


a 


-23.66 


10.37 


10.34 


2.89 


F 


27 


45 


aa 


2.3 


NGC 5457 (M101) 


SABcd 


6.96 


a 


-23.26 


10.23 


29.19 


5.19 


E 


37 


18 


dd 


5 



References. 
et al.l (120081). 



Dis tances: a = [ Jacobs et al. | <[2009b; b = [ Herrmann et al.l J2008j) ; c = jTulTy] | |1988ll; d = |Tully et all <|2009[|: e = |Tully| 
R25: A = Boroson (1981); B = Uarignan (1985); (J = Elmcgrccn & Elmegrccn (1984); U = Elmcgrccn & Elmcgrccn (1985); E = 
Koopmann ct al. (2006); F = Simicn &; dc Vaucoulcurs (1986); G = van Zee ct al. (1998). j/PA: aa = HypcrLcda (Faturci ct al. 2003); bb = 
Fiosma et al.l I u9au1); cc = Kamphuis & Briggs (1992j ra(l = [K crmicutt & liarnctt (199 m; ee =|van Zee fe Fjryant ( 19 9WT 
FSresolin ct al. (2009b); i = |Bresolin et al.|ftjUU5) ; 4 = |Garnetf et a!7| | |iy97| ; 5~ 



Breso l m eCa T.n|20U9a): 2 
et al.Hiy85| ;T - 



- |J_ 

Iran Zee eraTpfTggSft 

Note. — Col.(l): Galaxy name. Col. (2): Morphological type, from NED. Col. (3-4): Distance, in Mpc, and reference 
magnitude down to 20 mag/arcsec 2 isophote, from 2MASS. Col. (6): 8— lOOO^rn luminosity, determined using the formula in Sanders & Mirabel 



= |Kenmc"utt et al.| 



±1 11 r egio ns: 1 = 
p003[ ); 6 = |McCall| 



Col. (5): iv-band 



l |l996| and infrared fluxes from NED. Col. (7): B-band optical radius from 25 mag/arcsec 2 isophote, in kpc, from HyperLeda 1 PatureI*^t^I^2003*f 
except for M101, which is from the RC3). Col. (8-9): B-band stellar exponential disk scale length, in kpc, and reference. Col. (1U-12J: Inclination, 
in degrees; position angle of the galaxy line of nodes, in degrees east of north; and reference for these. Col. (13): Reference for emission-line data. 



Table 3 

Gradient Fits 



Galaxy 


N reg 


A(dex/kpc) 


A(dcx/R 25 ) 


Intercept 


(1) 


(2) 


(3) 






(4) 




(5) 




Control Sample 


NGC 300 


27 


-0.1246 




-0 


664 




8.95 


NGC 628 


25 


-0.0410 




-0 


536 




9.10 


NGC 925 


14 


-0.0471 




-0 


568 




8.83 


NGC 1232 


22 


-0.0310 




-0 


451 




9.10 


NGC 2403 


25 


-0.0750 




-0 


689 




8.93 


NGC 2805 


15 


-0.0237 




-0 


342 




8.93 


NGC 2903 


17 


-0.0239 




-0 


357 




9.08 


NGC 2997 


14 


-0.0642 




-0 


741 




9.30 


NGC 3184 


17 


-0.0613 




-0 


590 




9.35 


NGC 5236 


38 


-0.0413 




-0 


427 




9.29 


NGC 5457 


25 


-0.0310 




-0 


899 




9.13 


Interacting Galaxies Sample 


UGC 12914 


20 


-0.0071± 


0.0033 


-0 


104± 


0.048 


9.13± 0.02 


Arp 256 NED01 


13 


-0.0142± 


0.0022 


-0 


193± 


0.031 


8.95± 0.01 


Arp 256 NED02 


25 


-0.0162± 


0.0014 


-0 


271± 


0.023 


9.05± 0.01 


UGC 312 


13 


-0.0114± 


0.0070 


-0 


133± 


0.081 


8.84± 0.05 


UGC 813 


20 


-0.0210± 


0.0040 


-0 


233± 


0.045 


9.00± 0.02 


UGC 816 


21 


-0.0189± 


0.0023 


-0 


252± 


0.030 


9.05± 0.02 


NGC 2207 


43 


-0.0124± 


0.0026 


-0 


273± 


0.057 


9.08± 0.03 


IC 2163 


19 


-0.0165± 


0.0036 


-0 


162± 


0.035 


9.17± 0.02 


Arp 248 NED01 


8 


-0.0368± 


0.0040 


-0 


512± 


0.056 


8.85± 0.02 


Arp 248 NED02 


11 


-0.0144± 


0.0020 


-0 


190± 


0.026 


9.03± 0.01 


NGC 3994 


25 


-0.0224± 


0.0072 


-0 


139± 


0.045 


9.08± 0.02 


NGC 3995 


41 


-0.0175± 


0.0030 


-0 


339± 


0.059 


8.86± 0.03 


NGC 7469 


12 


-0.0146± 


0.0089 


-0 


179± 


0.109 


9.09± 0.05 


IC 5283 


12 


-0.0249± 


0.0060 


-0 


246± 


0.059 


9.14± 0.01 


UGC 12545 


14 


-0.0323± 


0.0142 


-0 


381± 


0.168 


8.87± 0.08 


UGC 12546 


20 


-0.0395± 


0.0082 


-0 


421± 


0.087 


9.10± 0.03 



Note. — Col.(l): Galaxy name. Col. (2): Number of H II regions 
used in fit. Col. (3-4): Radial oxygen abundance gradient, in dex/kpc and 
dex/R25- Col. (5). Intercept of oxygen abundance vs. radius fit, expressed 
as 12+log(0/H). 



